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Abstract 

This paper treats nonrelativistic matter and a scalar field cj) with a mono- 
tonically decreasing potential minimally coupled to gravity in flat Friedmann- 
Lemaitre-Robertson-Walker cosmology. The field equations are reformulated 
as a three-dimensional dynamical system on an extended compact state space, 
complemented with cosmographic diagrams. A dynamical systems analy¬ 
sis provides global dynamical results describing possible asymptotic behav¬ 
ior. It is shown that one should impose global and asymptotic bounds on 
A = —V~^dV/d([> to obtain viable cosmological models that continuously 
deform AGDM cosmology. In particular we introduce a regularized inverse 
power-law potential as a simple specific example. 


1 Introduction 


Recent observations, such as those of the cosmic background radiation, show that 
A cold dark matter (ACDM) cosmology is remarkably consistent with observations, 
although some tensions remain, see Hiq and references therein. On the other hand, 
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the unknown origins and nature of dark matter and dark energy remain mysteries 
which generate a steady flow of models and theories. Nonetheless, current obser¬ 
vations should be taken seriously, and viable cosmological models should therefore 
presumably only deviate from ACDM cosmology marginally. Moreover, although 
successful, the ACDM model still needs to be observationally tested, which suggests 
that it should be continuously deformed. In addition, for a variety of reasons, it is of 
interest to consider a dynamical dark energy content described by a held theoretical 
model instead of a cosmological constant A or a dark energy huid. This leads to 
the question: What phenomenological conditions must be imposed on an effective 
classical held description of dark energy to continuously deform ACDM cosmology 
for an open set of observationally viable solutions? 

Although observational compatibility arguably comes hrst, there also exist other 
issues that are of interest such as various hne-tuning problems and if there exist ob¬ 
servationally viable models that emerge from, or at least have ties to, a fundamental 
theory. The rather modest purpose of this paper, however, is to 

(i) describe general features and determine asymptotic behavior of hat Friedmann- 
Lemaitre-Robertson-Walker (FLRW) models with a perfect huid and a scalar 
held 0, minimally coupled to gravity, with a positive monotonic potential V (0), 

(ii) shed some light on the above “observationally viable held theoretical ACDM- 
deformation issue” by considering some simple models of the above type with 
the equation of state for the perfect huid specialized to dust, representing 
nonrelativistic matter (mainly baryons and nonbaryonic dark matter). 

Since the literature about minimally coupled scalar helds and a perfect huid in a 
hat FLRW spacetime geometry is vast, we will start with only a few general sample 
references: the inhuential papers mmn and the recent review |H], and references 
therein. 

It is of interest for discussions and future purposes to hrst consider a perfect huid 
with a general barotropic equation of state pm = (7m — I) Pm, where pm and Pm > 0 
are the pressure and energy density, respectively, and 7m = 1m{,Pm), together with a 
scalar held 0, minimally coupled to gravity, with a self-interaction potential V{(j)). 
Thus, the Einstein and matter held equations for the present hat FLRW models can 
be written as (see, e.g., HlEli) 


— 2 *;^^ + ^(0) + Pm, 

(la) 

H = — ^0 -|- 'ymPm^ , 

(lb) 

0 = 0 + 3H^ + V^, 

(Ic) 

pm 

(Id) 


where = dV/d(f) and where overdots denote derivatives with respect to cosmic 
time t. Units are such that SirG = 1 = c, where G is the gravitational constant 
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and c the speed of light in vacnnm. It is assnmed thronghont the paper that the 
Universe is expanding, i.e., H > 0. 


Eqnation (Ic) can be henristically interpreted as an eqnation for a particle of nnit 
mass with a one-dimensional coordinate 0 , moving in a potential V ( 0 ) with a friction 
force —3H4>. Eqnation (la) shows that H can be expressed in terms of 0, <p, and 
Pm- Introdncing 0 as an independent variable therefore leads to a three-dimensional 
dynamical system for ((),cl),pm- Alternatively one can solve Eq. (Id) to obtain pm = 
Pm{ci) and consider the variables 0 , 0 , a as the dependent variables, where a obeys 
the eqnation d = aH, where H can be expressed in terms of 0, 0, a by means of (la). 

It is of some interest to introdnce effective eqnation of state parameters, dehned by 


P<p + P<p ^ 

P4> |02 + p^( 0 )’ 


Ttot — 


Ptot Ptot 
Ptot 


or, alternatively, w* = = 7 * — 1 (* = 0 , tot), where 


P4> = + V"(0), 

Ptot Pcf) “1“ Pm 7 


p<t> = 

Ptot P^ “1“ Pm' 


( 2 ) 


(3a) 

(3b) 


Adding a cosmological constant to a perfect flnid with a barotropic eqnation of state 
can be regarded as a problem with a single flnid with a changed barotropic eqnation 
of state. Replacing a cosmological constant in a general barotropic flnid description 
with a scalar held therefore adds the dynamical degrees of freedom, 0 and 0 , to pm 
(or a), which implies a complication, since this leads to a larger set of solntions 
with different histories. In the context of ACDM deformations this problem is 
somewhat alleviated by the following assnmption: To obtain an evolntionary history 
snch that primordial nncleosynthesis gives observationally compatible light nnclear 
abnndances, we will assnme that there exists a mechanism snch as inhation that sets 
initial conditions for the present models at high redshifts so that the matter content 
greatly dominates over the scalar held (dark energy) content in the early Universe. 


For simplicity we will from now on consider a linear eqnation of state with a constant 
7 m restricted to the range 0 < 7 m < 2 , which avoids bifnrcations at 7 m = 0 , 
corresponding to a cosmological constant, and 7 m = 2 , which yields a stih perfect 
hnid with the speed of sonnd eqnal to the speed of light; 7 m = 1 corresponds to 
a hnid withont pressnre, which will be the focns when describing and discnssing 
ACDM deformations, while 7 m = 4/3 corresponds to a radiation hnid. It follows 
from Eq. (Id) that 


pm — Poifl/ ®o) 


- 37 „ 


(4) 


Fnrthermore, most of the scalar held potentials that have been considered as possible 
alternatives to a cosmological constant as dark energy have strictly monotonic scalar 
held potentials, as exemplihed by, e.g., the popnlar inverse power-law potential. 
Apart from the model with U = A = const., we will therefore assnme that V is 
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defined, differentiable, positive, and strictly monotone for (j) G ((/>_, 0+). Withont 
loss of generality, we assnme that V is monotonically decreasing. We also assnme 
that the potential is snch that 0+ = oo. If 0_ is hnite or — cx) depends on how fast 
1/(0) increases when 0 —>■ 0_, e.g., for an inverse power-law potential, V oc 0“", 
0_ = 0, while an exponential potential, V oc exp(—A0), leads to 0_ = —oo. 

A constant or monotonically strictly decreasing potential makes it convenient to nse 
a particnlar type of dynamical systems formnlation that brings 

m = -y (8) 

into focns, where A is zero for a constant potential and positive for 0 G (0_,oo) 
when V (0) is strictly monotonically decreasing. Within this context it is natnral to 
regard A as more “fnndamental” than V itself, where V is obtained from A(0) via 

V = Vbexp f— f X{^)d4) 

V JlPo 

In the special case where A and hm(^^ 0 _ A are bonnded it follows that the potential 
can be bonnded by an exponential and therefore 0_ = —oo. If A is nnbonnded, as 
in the case of an inverse power-law potential, 0 _ is hnite. 

Finally we note that the system ([^ has nnbonnded variables and right-hand sides 
that blow np. It is therefore not snitable for a global or asymptotic analysis of the 
solntion space and its properties. Below we will therefore make a change of variables 
in order to obtain a dynamical system on a compact state space that allows these 
issnes to be addressed, as well as giving global illustrative pictures of the entire 
solution spaces of the models under consideration]^ 

The outline of the paper is as follows. In the next section we (i) formulate the 
dynamical systems approach that is used throughout the paper to deal with a huid 
and a minimally coupled scalar held; (ii) describe the state space features; (iii) point 
out some global properties and determine asymptotic behavior for models with a 
perfect huid and a scalar held such that 0 < A(0) < cxd, 0 G (0_, oo), and 


( 6 ) 


lim A = A_|_, lim A = A_, (7) 

where A+ is hnite and A_ is hnite or oo. 

In Sec. we (i) situate ACDM cosmology, which occurs for A = 0, in the three- 
dimensional state space, which leads to the concept of an attracting separatrix 
surface rather than individual “attractor solutions”; (ii) consider the A = const, 
(i.e., an exponential potential) and dust models in the present context and give a 
specihc example that shows that some of these models are completely solvable; (iii) 

^For a few examples of work dealing with dynamical systems methods in cosmology, see dill 
[iniiiiiiiaini Halil], and references therein. 
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discuss the symmetries of the above models and how the associated structures might 
be somewhat misleading for the general picture. 

In Sec.j^we turn from the previous “frozen A” models to some “dynamical A” models. 
We consider (i) the inverse power-law potentials, characterized by A = a/0, in a 
global dynamical systems setting, tying local results to the global features discussed 
in Sec. in particular we point out that observational viability requires hne-tuning 
of a to small values, since A —)■ oo when 0 —)■ 0 produces a “memory” that affects the 
entire evolution of, e.g., “tracker” (attractor) solutions. For this reason we therefore 
(ii) consider the arguably simplest possible “A-regularization” of the inverse power- 
law potential, namely, A = aj which, in contrast to the inverse power- 
law case, gives rise to continuous deformations of the attracting ACDM separatrix 
surface, and which for C > a yields much less stringent observational compatibility 
conditions on a. Finally, the section is concluded with (iii) a general discussion 
about observational viability conditions that A should satisfy for any potential that 
might arise from some fundamental theory. 

Throughout we use dynamical systems pictures which clearly display general key 
features, and, in order to compare with ACDM cosmology, these pictures are com¬ 
plemented by diagrams that indicate the evolutionary history of physically important 
quantities. Finally, we outline how it is possible to treat two fluids, e.g. dust and 
radiation, and a scalar held in the Appendix. 


2 Dynamical systems approach 

2.1 Dynamical systems formulation 

The presently used dynamical systems formulation is based on the following depen¬ 
dent variables, 

( 8 ) 

These variables are by no means suitable for all scalar held potentials, e.g., for 
potentials with a zero minimum such as monomial potentials it is advisable to replace 
the scalar held variable Z with an FT-based variable that takes into account a varying 
averaged oscillatory time scale at late times, as done in [15]. Nevertheless, the above 
variables are useful for positive monotonic potentials. An optimal choice of the 
variable Z (0) depends on the potential that is studied. However, to obtain a suitable 
global dynamical systems formulation Z should always be chosen to be a globally 
invertible monotone function in 0, dehned on a bounded interval Z G (Z_,Z+), 
Z± = Z{(j)±); without loss of generality, we choose Z to be monotonically increasing 
in 0 so that dZ/dcj) > 0 for Z G (Z_, Z+). 

Apart from the dependent variables we also need to choose a new time variable. In 
order to obtain a regular dynamical system this choice depends on if A is bounded 
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or not. In the case A is bounded, as exemplified by, e.g., exponentially bounded 
potentials, it is convenient to use r, defined by dt = H~^dT, where r = ln(a/ao) 
sometimes is referred to as N, the number of e-folds from a reference time to [oq = 
a(to)], which for the present time leads to r = — ln(l -|- z), where z is the redshift. 
Thus, the field equations can be written as the following coupled system: 


-{2-q)x+ Y^A(Z) Gy, 

(9a) 

3 [2x‘^ -7m(l -G^)] flm, 

(9b) 

Ve—x 

(9c) 


where ' denotes derivatives with respect to r. With some slight abuse of notation 
we have written A(Z) = A(0(Z)) = —V^jV. Furthermore, the Gauss constraint (la) 
is used to globally express the Hubble-normalized scalar field potential energy, Gy, 
in terms of the state space variables: 

= ( 10 ) 

while the deceleration parameter g, defined via H' = is given by 


q = -I+ ?,x^ P 


( 11 ) 


The above system is suitable when A is bounded, but not when lim(^^<^_ A = oo, 
as exemplified by, e.g., inverse power-law potentials. In order to obtain a regular 
dynamical system when lim^^^_ A = cxd we therefore choose a new time variable r, 
defined by dr = g[Z)df, where g{Z) is a suitable positive bounded function of Z 
such that lirciz^z- = 0) 

lim gx{Z) = gx_ = const. < oo, where gx{Z) := g{Z)\{Z), (12) 


and \\mz^z+ g{,Z) = 1, where, e.g., g{Z) = 1/(1 -|- X{Z)) is a suitable choice for 
inverse power-law potentials. This results in the system 


d T" / 

- = -g{Z){2 - q)x + sj-gx{Z)Dy, (13a) 

= ^g{Z) [2x^ - 7m(l - G,^)] G^, (13b) 

f^y6g(Z)|.. (13c) 

Note that @ is obtained by setting g{Z) = 1, hence r = r, and so the above system 
formally covers both cases. 
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2.2 State space structures 


The above assumptions lead to the following domain for the state vector (x, flm, 
for the interior state space S for models with a perfect fluid and a scalar held: 

rim > 0, riy = 1 - - rim > 0, Z_ < Z < Z+. (14) 

Furthermore, it follows from ([^, and the following auxiliary equation for riy, 

n'y = 2(^l + q- yf A(Z) a; j (15) 

that it is possible to extend the state space and include the following boundaries: 


The boundary rim = 0 (i.e. pm 
subset [see Eq. (9b)], and thus 


0) yields the (pure) scalar field boundary 


+ riy — 1, 


(16) 


with an interior state space S<^ given by riy = 1 — > D, Z_ < Z < Z^. 

• The boundary riy = 1 — — rim = 0 corresponds to a model with a perfect 

huid and a massless scalar held; we refer to this subset as the massless scalar 
field boundary subset. Note that this boundary also contains the pure massless 
scalar held boundary subset rim = 0, x = ±1, and the perfect fluid boundary 
subset rim = 1 (ri0 = 0, i.e., riy = 0 = x), which is described by a line of hxed 
points, FL^, with constant Z. 


• In addition we assume that the potential and Z are such that the state space 
S can be extended to not only include the boundaries rim = 0, riy = 0, but 
also the boundaries Z = Z_ and Z = Z^, which furthermore are assumed to 
constitute invariant boundary subsets. 


Thus the extended compact state spaces S and S</, are characterized by 

fim > 0, fly = 1 - X^ - flm > 0, Z_ < Z < Z+, (17a) 

= 1 - a;2 > 0, Z_<Z<Z+, (17b) 

respectively, and the dynamical systems on S and are of diherentiability class C^ 
or higher. We therefore assume that the potential is such that it is possible to hnd 
a scalar held variable Z and, if needed, a function g{Z) so that 


(i) g\{Z) and g{Z)dZ/dfi are bounded and diherentiable on the extended interval 

ze [z_,z+], 

(ii ) g{Z)dZ/d(f)\z=z± = 0, i.e., Z = Z_ and Z = Z+ are invariant boundary 
subsets. 
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Although the above restrictive assumptions require that the potential 1^(0) is of 
differentiability class they cover a vast class of potentials V{(()). Moreover, 
the extended state space treatment of the admissible class of potentials can be 
used to provide bounds on asymptotic behavior of potentials that do not admit Z- 
extended state spaces, and thus they provide a natural starting point for very general 
situations. It should also be noted that for many qualitative dynamical aspects, as 
will be seen, there is no need to state any more information about Z{(j)) or 
although quantitative treatments of course need specific functions and thus we will 
give specific examples when dealing with continuous ACDM deformations. 

To physically compare solutions that are obtained in the dynamical systems ap¬ 
proach with ACDM cosmology, we will complement the state space description with 
diagrams involving physically important quantities. Apart from the physically im¬ 
portant state space variable we have given the deceleration parameter q in terms 
of state space variables in Eq. 0 . and below we will also express H in state space 
variables. In addition it is of interest to give the jerk parameter j [MlEI , since even 
though the observational feasibility of j is questionable, it is still useful in order to 
theoretically compare models with ACDM cosmology for which j = 1; for constant 
7m the jerk parameter is given by 

j = l + Qx^ - ^7m(l - lm)Om “ ^'/Q\{Z) Dy X. (18) 

Next we turn to structures that are helpful when it comes to determining global and 
asymptotic state space features. 


2.3 Nonlocal features 


The monotonicity of H 

The Hubble parameter H, or more conveniently can be expressed in the state 
space variables in S and and is given by 


2 _ V{cj>{Z)) _ V{cj>{Z)) 


= 


SQv 3(l-x2-Dm)’ 


(19) 


as follows from the definition of Dy. This has far-reaching consequences for the 
global dynamics due to 

7 (H^y = -2{l + q)H^, 


H' = -{l + q)H 
since Eqs. (10) and ([II|) lead to 


( 20 ) 


1 + <? = + ^7mDm, 2 - g = 3Dy + ^(2 - 7m)Dm, (21) 


which together with 0 < 7 m < 2 yield 

- 1 < g < 2. 


( 22 ) 
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As will be seen, the above inequalities for q, induced by energy and causality con¬ 
ditions, play a central role for the dynamics, which strongly suggests that the held 
equations should be expressed in terms of q, as done in (|^. 


The above relations imply that q = 2 only occurs on the pure massless scalar held 
subsets X = ±1, while q = —1 requires Dm = x = 0, fly = 1. Since 1 -|- g > 0 
on S, is a monotonically decreasing function on S. This is also true for S^, 

as can be seen as follows. On we have the inequalities —l<g< 2 . As a 
consequence is a monotonically decreasing function also on this subset, except 
when 1 -|- g = 0. However, when 1 -|- g = 0 is not an invariant subset, then 1 -|- g = 0 
just corresponds to an inhection point in the graph of and is thereby still 
monotonically decreasing]^ It is only when 1 -|- g = 0 is an invariant subset that 
stops decreasing and for which 1 -|- g = 0 can be an asymptotic state. For this to 
happen requires that x = 0, since 1 -|- g = 3x^ on S^, also means that x = 0 must be 

an invariant subset. Since x' = when fl^ = 0 = x, fly = 1, this leads to the 

condition that X{Z) = 0, which by our assumptions about monotonically decreasing 
potentials can only happen on the boundaries. Hence is also monotonically 
decreasing on S<^. This implies that there can be no hxed points or recurring orbits 
(solution trajectories) in S and S^, and due to (19) all solution trajectories in S and 
originate and end at Z± and fly = 0, which shows that it is essential to extend 
the state space to include these boundaries in order to fully describe the dynamics 
of the present models . 


The Z± boundaries 

Let us begin by considering the Z± boundaries for the case where \imz^z± A = A± 
is hnite. Dropping the subscript ± on A± leads to the system 

x'=-(2 - g)x-f ^^Afly, (23a) 

= 3 [2x2 - 7^(1 - D^)] D™, (23b) 

on Z± where we recall that fly = 1 — x^ — fl^ and g = — 1 -|- 3x2 _|_ As will 

be elaborated on in the next section, these equations describe the reduced equations 
of an exponential potential V oc exp(—A0) when A 7 ^ 0 and those of a constant 
potential when A = 0. The above system has an extended state space Sred given by 

f^m > 0, Dy = 1 - x2 - fl^ > 0, (24) 

and admits a number of hxed points given in Table 

We have here introduced a notation for the hxed points where the kernel M stands 
for a massless scalar held state; FL stands for a Friedmann-Lemaitre perfect huid 
state, while dS stands for a de Sitter state (for interpretation of various types of 


^For an elaboration of this in the case of monomial potentials, see m- 
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FL 

X = 0, flrti = 1 

q = (37m - 2)/2 

- 

EM 

X = y/3/2'ymX~\ f2m = 1 - 37mA“/ A > y/S'y^ 

q = (37m - 2)/2 

7</> ~ 7m 

M± 

X = ±1, Qm = 0 

q = 2 

7^ = 2 

PL 

X = A/ y/E, fim = 0; A < y/E 

g = -1 + AV2 

70 = AV3 

dS 

X = flrn = 0; A = 0 

q = -1 

70 = 0 


Table 1: Table depicting the type of fixed points that can occur for a FLRW model 
with a monotonically decreasing scalar held potential with a hnite A and a perfect 
huid with a linear equation of state parameter 0 < 7 m < 2. 


de Sitter states in scalar held cosmology, see [IS]); the kernel PL stands for power 
law (inhation, when A < \/2, since this leads to an accelerating state with g < 0), 
while EM stands for a scaling solution with an exponential potential and matter in 
the form of a perfect huid with a linear equation of state. In addition a superscript 
describes the value of x, while in the full three-dimensional treatment we also add 
a subscript that describes the value for Z. 

A one-parameter set of orbits in Sred originates from each of the sources M"*" and 
M“ when A < y/E, but when A > \/6 there are no solutions that originate from 
M"*" into Sred, while one solution originates from FL. Toward the future PL is a 
sink when A < y/S'jrn while EM is future stable when A > y/3'ym- Thus bifurcations 
take place when A = a/6, which is when PL leaves the physical state space via M"*", 
and when A = -/37m, which is when EM enters the state space via PL. In addition 
there is a bifurcation when A = 0, which is when PL becomes the future sink dS. 
The solution space structures for representative examples of these different cases are 
given in Fig. 

The case when A —)■ oo when / —yields the following system on Z_: 



df 


(25a) 

(25b) 


where we recall that g\_ = hm^^^_ (g(Z)A(Z)) > 0. It follows that flm = const, 
while the subset Qy = I — x^ — fim = 0 forms a line of hxed points M^, and that 
X is monotonically increasing in Sred- Thus M^ with a; < 0 is the source while M* 
with X > 0 is the sink of all orbits in Sred; hnally, note that M* with x = 0 is the 
hxed point FL. 


The Dy = 0 boundary 

We hrst note that the equation for x on fly = 0 is invariant under the transformation 
X —>• —X, which enables one to use flstis = x^ as a variable instead of x, and as a 
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(c) Solution space for A = 2 and 7 ™ = 1. 


(d) Solution space for A = 3 and 7 m = 1. 


Figure 1: Representative solution space structures for Sred- 


result the dynamics of the massless scalar held explicitly takes the same form as 
that of a stiff huid, p = p- This leads to the equations 






Z' = = ±V6 ^ \/l-VLr^ 

dcp dcp 


(26a) 

(26b) 


where we have used the Gauss constraint to write Gstifr = 1 — Dm and x = 

Since the equation for Dm decouples from that of Z, it follows that all orbits (solution 
trajectories) in the two-dimensional state space look the same as the orbits of the 
one-dimensional problem associated with Dm, when projected onto Dm- As follows 
from (26a), Dm monotonically increases from 0 to 1. Thus the past asymptotic state 


of all orbits on Dy = 0 resides on the pure massless scalar held subset Dm = 0, 


x^ = 1 and thereby q = 2] thus either x = -|-1, which from (26b) leads to Z = Z_ 


and thus a hxed point M^_, since dZ/dcj) > 0, or x = —1, which yields Z = Z^ 
and a hxed point It also follows that the future asymptotic state of all orbits 
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on lly = 0 resides on 0^ = 1, x = 0, Z = const., which thereby form a line of 
hxed points FL^. To determine what happens with Z toward the future on fly = 0 
we note that it is possible to solve the equation for Qm explicitly, but the solution 
can also be obtained by using the fact that pstifi oc exp(—6r), pm oc exp(—3r), 
= pstiff + Pm, which gives 


TT = (f^stifroexp(-6r) + exp(-37mr))^'^^ , 

Ho 


(27a) 


X = ±fl 


1/2 

stifFO 


exp(-3r)/(i7/ifo), 


(27b) 

= fimoexp(-37^r)/(i7/i7o)^- (27c) 

To obtain the solution for Z requires a specihcation of Z. However, it is easier to 
obtain Z by integrating the equation for 0, 

0 ' = x/ 6 a: = ±(6f2stifro)^'^^ exp(-3r)/(i7/i7o), (28) 

and then inserting the result in the chosen dehnition of Z. As follows from the 
above equation, 0 —)■ const, when r —)■ +cxd, and therefore also Z —)■ const, in this 
limit; i.e., each point on the line FL^ attracts one solution with positive and one 
with negative initial x on the fly = 0 boundary. For a representative depiction of 


the fly = 0 boundary, see Fig. 3(a) 


Role of the fly = 0 boundary 

To study the neighborhood of fly = 0 we note that 

^y^^ybv=o = 2 + g — \j~^X{Z) , (29) 

and hence 

T A(Z^)) , fly'fl'ylFL^ = 37m. (30) 

As a consequence it follows that is a source for orbits in S and S<^ as is M)^_ 
when A(Z_) < Vg, but if A(Z_) > \/6, then no orbits enter S or S<^ from M^ j^a 
one-parameter set of orbits, one from each point on FL^, enters S, where FL^ acts 
as a “transversal saddle line,” forming a two-dimensional invariant subset in S. By 
demanding initial data for which fl^ is close to 1, this implies the important feature 
that for all the models we consider such solutions can be approximated by this (at 
least initially) attracting invariant surface subset of codimension 1 with respect to 
the state space dimension. We will later contrast this feature with discussions about 
single “attracting” solutions that appear for various models in the literature. 

^Solutions originating from corresponds to “scalar field particles” that come from oo 

moving toward decreasing (j), while if A(Z_) < \/6 the potential is shallow enough so that scalar 
field particles also can come from —oo moving toward increasing b- Solutions originating from 
FL^ correspond to scalar field particles that initially lie still at some initial value b and then roll 
down the potential F(b). 
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Monotonicity of negative x 


Further restrictions on the global dynamics of the present models come from eqs. (9a) 


and (9c). It follows from (9a) that x is monotonically increasing when — 1 < a: < 0, 


since then 2 — g > 0. This in turn implies that if an orbit resides or partly resides 
in the region —1 < a: < 0, then the orbit originates from x = —1, which due to (9c) 


and dZld(f) > 0 implies that the past asymptotic state must be Furthermore, 

the monotonicity of x when — 1 < x < 0 implies that all orbits in S and must 
end at X > 0 . 


The equation for x can be written on the following form: 






X 


FT j 


(31) 


which shows that x is monotonically decreasing if 1 > x > X{Z )/ a/6. In combination 
with the above result for x this implies that all orbits in S and must end at 
0 < X < A(Z)/a/6. However, in order for an orbit to end at x = 0 requires that 
X = 0 be an invariant subset. Since 


X |a:=0 — 



A(Z)(1-SJ„), 


(32) 


and since flm = 1 is the FL^ subset, it follows that this can only be the case 
if A(Z) = 0, which can only happen on the Z± boundaries for a monotonically 
strictly decreasing potential. The above equation also shows that x is monotonically 
increasing when x = 0 when < Z < Furthermore, since dZ/dcp > 0 it follows 


from (9c) that Z is monotonically increasing (decreasing) in r, or r, when x > 0 
(x < 0 ). 

If X is future asymptotically positive it follows that Z —)■ However, to exclude 

or show that any of the orbits from end up at dS 2 _ if A_ = 0 requires more 
information about the potential. If lim 0 ^^_ V{(j)) = const., then there are solutions 
that end at Z_, but if H —)■ cxd when 0 there are not. This can be understood 


by considering (Ic). A solution with negative x can be viewed as a particle moving 
in a direction where the potential is monotonically increasing. At the same time 
it is losing energy due to the friction force —SHcp until 0 = 0, but then 0 = —I/,, 
where is strictly negative. As a consequence the particle changes direction and 
0 and thereby also x become positive, and therefore Z —)■ Z+ toward the future. 


Thus all orbits in S and end at the invariant boundary subset Z_|_ when the 
potential is monotonically decreasing and H —?• oo when 0 —>■ 0 _, which we for 
simplicity from now on assume (apart from the case of a constant potential, which 
is treated separately, all our explicit examples obey these conditions). 
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The role of the Z_ subset 

The four dynamically qualitatively different cases on Z_ when 0 < A_ < cx) lead to 
different interior dynamics: 

(i) When A_ < v^3^, a two-parameter (one-parameter) set of orbits originates 
from the source M^_ into S (S<^), while a single orbit originates from PL^_ 
into S^. 

(ii) When y/S'jm < A_ < a/6, again a two-parameter (one-parameter) set of orbits 
originates from the source M^_ into S (S^), but in this case there is also a 
one-parameter set (a single orbit) that enters S (S^) from PL^_ as well as a 
single orbit from EM^_. 

(iii) When A_ > \/6, a single orbit enters S from EM^_. 

(iv) When A_ = 0, a two-parameter (one-parameter) set of orbits originates from 
the source M^_ into S (S,^), while a single orbit originates from dSz_ into S^. 

When A_ —)■ CX) at Z_, there is just one orbit that enters S from FL^_, the so-called 
tracker or attractor solution. The exclusion of orbits into S and S<^ from the fixed 
points M^_ with a; / 0 can be understood as follows: There are no solutions coming 
from Mf_ when a; < 0 since Z is decreasing toward the future when a; < 0. When 
a; > 0, it follows that fly is decreasing in the vicinity of Mf_ and hence there are 
no solutions coming from this part of M^_ into S or either. Hence M^_ with 
a; = 0, i.e., FL^_ is the only fixed point on M^_ which can give rise to a solution 
into S, and it does give rise to the tracker solution, which can be established by 
center manifold theory or by consideration of nearby solutions and continuity, since 
this fixed point acts as a kind of saddle in the full state space. 


The role of the Z^ subset 

The only solutions that originate from the Z^ subset into S (S<^) are those from the 
source from which a two-parameter (one-parameter) set of orbits originates 

into S (S,/,). 

For the present monotonically decreasing potentials, all orbits in S and end at 
as follows when 0 < A+ < cx: 

(i) When A+ < y/3^, a two-parameter (one-parameter) set of orbits ends at the 
global sink PL^^ from S (S^). 

(ii) When y/S'jrn < A+ < ^/6, a two-parameter set of orbits in S ends at the sink 
EM^^ while a one-parameter set of orbits in end at PL^^. 
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(iii) When A+ > ^/6, a two-parameter set of orbits in S ends at the sink 
while a one-parameter set of orbits in end at 

(iv) When A+ = 0, a two-parameter (one-parameter) set of orbits ends at the global 
sink dSz^ from S (S<^). 


An additional useful quantity 


Finally, consider 


which obeys 





(33) 

(34) 


as follows from (9b) and (15). 


When 



X(Z) X > 0, ^ is monotonically 


decreasing. For globally bounded A such that X{Z) < it follows that ^ is 

strictly monotonically decreasing and hence 0 and thereby t 0 toward the 
future; i.e., the future dynamics resides on the scalar field subset. Toward the past 
^ oo and hence fly —)■ 0 and as a consequence it follows that the past dynamics 
asymptotically resides on the fly = 0 subset, which of course is already implied by 
our previous discussion in a somewhat more complicated manner, which, however, 
includes more details and general situations. 


We will now turn to explicit examples of potentials and a fluid with a dust equation 
of state, 7m = 1- Throughout we will choose scalar field variables so that = 0 
and = 1, so that S is given by 

flm > 0, fly = 1 - - flm > 0, 0 < Z < 1. (35) 


It should be pointed out that this type of choice of scalar field variable is by no means 
an optimal choice for all scalar field potentials that can be dealt with by means 
of a scalar field variable Z(0), as illustrated by the treatment of the generalized 
Chaplygin gas in m, where —1 < Z < 1, but choices such that 0 < Z < 1 
are “good enough” for our present purposes. Since Z_ = 0 and Z+ = 1 we will 
introduce 0 and 1 as subscripts instead of Z± to distinguish fixed points on the two 
boundaries (in addition we have the subscript Z for the line of hxed points FL^ and, 
for a constant potential, the line of hxed points dS^). 


3 Frozen ACDM dynamics 

Before considering the constant and exponential potentials, let us hrst recall some 
facts about ACDM cosmology. In their simplest form, the ACDM models are so¬ 
lutions of Einstein’s equations with (i) a spatially isotropic and homogeneous hat 
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FLRW geometry and (ii) a perfect fluid that has negligible pressure, i.e., dust, repre¬ 
senting nonrelativistic matter, and a cosmological constant A. It is useful to describe 
solutions by the following quantities: 


H = a/a. 





= 


Pm 

3^’ 

A 

3^' 



-q' + q + 2q^, 


(36a) 

(36b) 

(36c) 

(36d) 

(36e) 


Here the cosmographic (or cosmokinetic) quantities a(t), H(t), q(t), j(t) are the cos¬ 
mological scale factor and the time-dependent Hubble, deceleration, and jerk pa¬ 
rameters, respectively [lain], which take the following form for the ACDM models: 


-fT = V^moexp(-3r) IIao, 

-no 


Qrno exp(-3r) 


3 3 

^ 2 2 exp(—3r) -|- Hao J ’ 

j = 1, 

^ _ ^mo exp(-3r) 

fImoexp(-3r)-FII ao’ 


Ha — 


n 


AO 


fImoexp(-3r) IIao’ 


(37a) 

(37b) 

(37c) 

(37d) 

(37e) 


where + IIa = 1 and where Hq, flmO) and Hao are the values of H, and Ha 
at the reference time to and hence r = 0|^ Throughout we will use the following 
present values: 

f^mo = 0.3, Hao = 0.7. (38) 


It follows from ( |37b ) that go = —0.55. 

Figure [^depicts the evolution of H/Hq, q, and Dm in terms of the redshift 2 ; for the 
ACDM models with the above present values. As can be seen. Dm monotonically 
decreases from 1 (corresponding to the big bang limit r —>■ — cxd) to zero (corre¬ 
sponding to the inhnite future limit r —)■ -|-oo, described by a de Sitter state with 


“^Due to noisy data, the current observational feasibility of the jerk parameter j is questionable. 
Nevertheless, the fact that j = 1 for the ACDM models makes it useful as a theoretical indicator of 
deviations from ACDM cosmology, and it is for this reason we give j. Moreover, as will be seen, the 
jerk parameter seems to be the parameter that is most sensitive when it comes to deviations from 
ACDM cosmology; it is unfortunate that it is also the most difficult quantity to observationally 
measure of the presently discussed ones. 
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-1 012345678 

z 

(a) i?(z)-diagram for the ACDM models. 




(b) g(z)-diagram for the ACDM models. 


(c) nm( 0 )-diagram for the ACDM models. 


Figure 2: Depiction of the evolutionary history of ACDM cosmology for Q^o = 0.3, 
f^Ao = 0.7 in terms of the redshift z. 


Da = 1). Alternatively this also follows by considering the differential equation for 
Qfn, which is given by (see, e.g., p. 62 in 0 ) 


= -3(1 - 


(39) 


3.1 ACDM dynamics as A = 0-CDM dynamics 

Setting Pm = 0 and V = A yields a model that can be interpreted as that of dust, 
a massless scalar held, and a cosmological constant A. The associated model does 
not solve any of the issues one might hnd unattractive with the ACDM model, 
but it does show some aspects that are important in a broader context. Since 
A = 0, one obtains a reduced coupled two-dimensional system for x and D^,, while 
the decoupled equation for Z can be integrated once the evolution for x and Qm 
is found. However, instead of just considering the essential reduced extended two- 
dimensional state space Sred, we need to consider the full three-dimensional extended 
state space S to illustrate how this model is connected with more general models for 
which A = X(Z). Note, however, that the decoupling of the equation for Z from the 
coupled system of x and Qm results in that all solutions with the same initial data of 
X and Qm yield the same curves when projected onto the x — D^-plane, irrespective 
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of the initial value of Z. To proceed requires that Z be specihed. There are many 
possible choices that lead to an analytic dynamical system on a compact state space 
that cover the domain = ±cx), but here, due to future purposes, we choose 


1 + exp{—X(j)) ’ 

which is monotonically increasing in (j) and where A is an arbitrary positive constant; 
the scalar held variable Z thereby has an extended range Z G [0,1]. This leads to 
the following dynamical system for the present models: 


(40) 


x' = —{2 — q)x, (41a) 

= 3 [2x2 - (1 - (41b) 

Z' = v^AZ(l - Z)x, (41c) 


where 2 — q = 3(1 — x^) — = 3(1 — x^ — Qm) + \^m- The reduced coupled 

system for x and Elm is invariant under the transformation x —)■ —x, which makes it 
possible to identify this problem with that of a stiff huid, characterized by flstis = 
a cosmological constant, and dust. In a similar manner as for the = 0 subset, 
this leads to 

-JT = (hlstiffo exp(-6r) + IIao + exp(-3r))^'^2 ^ (•42a) 

Do 

X = ±V^stifroexp(-3r)/(hf/hfo), (42b) 

Drn = Emo exp(-3r)/ {H/Hof. (42c) 


As for the fly = 
given by 


0 subset it is easiest to obtain Z by integrating the equation for 0, 
(j)' = \/6x = ±(6flstifro)^^^ exp(-3r)/(Tr/Tro), (43) 


and then inserting the solution into (40). 


We have previously dealt with the global dynamics of monotonically decreasing 
potentials. We here give a complete description for the constant potential. The past 
and future states of all orbits follow from the exact solution, or from the monotone 
function in combination with the dynamical structure at the boundaries and 
the local properties of the hxed points. The result is that a two-parameter set of 
orbits in S and a one-parameter set in originate from each of the hyperbolic 
sources Mq and M)," while a one-parameter set of solutions originates from the 
transversally hyperbolic line of hxed points FL^ into S; orbits in S and S<^ end at 
the transversally hyperbolic line of hxed points dS^, which thereby constitute the 
future attractor. The orbits that come from M(J" (M^) correspond to solutions with 
increasing (decreasing) 0 from 0 —)■ —00 (0 —>■ -|-oo), starting from a massless state 
[heuristically this can be understood from Eq. (0 where the friction force —3H(j) 
generates energy toward the past, making the energy content of A (and p) negligible 
compared to the kinetic energy of the scalar held], while solutions that originate 
from FL^ reside on the invariant subset x = 0. 
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Since the x = 0 subset corresponds to solutions with different constant 0, and 
thereby also a constant = A, these models are identical to the ACDM models. 
We therefore denote the x = 0 invariant subset surface as the ACDM surface or 
subset. Furthermore, this surface divides the state space into two disconnected 
domains with x > 0 and x < 0, respectively, and the ACDM subset thereby acts 
as a separatrix surface. Note that it is a necessary condition that the potential is a 
constant for x = 0 to be an invariant subset in S and thus no other potentials yield 
a ACDM model. 

Due to the regularity of the dynamical system and that each fixed point on FL^ 
acts as a transversal saddle point (exemplifying the previous general discussion con¬ 
cerning the Dy = 0 subset and the vicinity of Dy = 0 and FL^), and since dS^ 
is the future attractor, it follows that the ACDM separatrix surface x = 0 is an 
“attracting” surface or invariant subset (indeed, x^ is monotonically decreasing). 
As a consequence, demanding that Qm is close to 1 initially leads to solutions with 
almost ACDM behavior. The solution space for this set of models is depicted in 
Fig. I where we have chosen A = \/3/2. Note that the structure on the Dy = 0 
subset is qualitatively the same for all the potentials we consider, as discussed in 
the previous section. 

The present models in effect do not really have a dynamical dark energy since the 
dark energy content “freezes” toward an “unnatural” value A > 0 toward the future 
for models with x 7 ^ 0 , and they therefore do not solve any coincidence or energy 
scale problems one might have with the ACDM models. Nevertheless, since models 
with Elm close to 1 as an initial condition basically have the same evolution as the 
ACDM models, these models with their attracting ACDM separatrix surface subset, 
in an evolutionary history sense, are the best a scalar field and dust model can 
accomplish when it comes to mimicking ACDM cosmology, and they therefore, from 
a purely time development perspective, set the standards for any truly dynamical 
dark energy model. Next we turn to the well-known case of an exponential potential 
(see, e.g., mm), but in the present three-dimensional context. 

3.2 Constant A ACDM dynamics 

An exponential potential 


V = Voexp(—A0), A > 0, 


(44) 


and a dust matter equation of state results in the dynamical system 



(45a) 

(45b) 

(45c) 


= 3 [2x2 - (1 - nj] 
Z' = V6AZ(1 - Z)x, 
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(a) The massless scalar field boundary fly = (b) The scalar field boundary 11^ = 0. 

0 . 


'YZ 


MT 



tential and dust. 


Figure 3: Depiction of the solution space of a fluid without pressure and a scalar 
held with a constant potential D(0) = A > 0, where Z = (1 + exp(—-^3/20))“^, 
illustrating the x = 0 ACDM separatrix surface subset and its attracting nature. 


where we have used the dehnition (40) for Z (also, recall that 2—q = 3(1—x^) — 

Once again there is a decoupling between the scalar held variable Z and a coupled 
system for x and Qm on a reduced state space Sred, parametrized by A. As for a 
constant potential, it therefore follows that solutions with the same initial x and 
Dm, but with diherent initial Z, have the same trajectories when projected onto the 
X — Dm-plane. 


The qualitative structure of the solution space of S is entirely determined by the 
general monotonic features and the structures on the boundary subsets of S and their 
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neighborhoods described earlier. Thus there are two-parameter sets of orbits in S 
that originate from the sources and M]" when A < a/6, but only from when 
A > \/6, while a one-parameter set of solutions originates from the transversally 
hyperbolic line of fixed points FL^. Toward the future PLi is a sink when A < \/3 
(and the future attractor on S and S<^) while EMi is future stable when A > \/3 
(and the future attractor on S). As previously shown, it is a general feature that 
there is a one-parameter set of solutions that originate from the line FL^, forming 
an attracting invariant subset, but it is only when A < a/3 that this invariant subset 
forms an “attracting separatrix surface,” which in the limit A —)■ 0 leads to the 
previous ACDM surface. 

It is of historical interest to point out that it was the case with A > a/6 with 
EMi as the future attractor (although in the two-dimensional projected context) 
that was the original role model for concepts such as attractor (the hxed point 
EMi is the formal future attractor) and tracker solutions in scalar field cosmology 
[recall that 7 <^ = for general (see. Table [^; i.e., the effective scalar field 
equation of state “tracks” the matter equation of state]. As we will see, this is 
rather misleading. Moreover, from current observational constraints these models 
are no longer viable cosmological models since q = 1/2 for the future attractor 
EMi. Since models with small A are arguably more relevant for understanding more 
general and observationally competitive scalar held models, we focus on this class 
of models next. 


It is worthwhile to point out that since this does not seem to be generally known, 
there are models with a perfect huid and a scalar held with an exponential potential 
that admits simple explicit solutions, as shown by Uggla et al. in |20] (where many 
other solvable scalar held and modihed gravity models can be found as well). One 
such example is for dust when A = \/3/2. Using the methods in [20], we see that x 
and Qrn can be found in terms of cosmological time t\ 


X = 


2t^ - 


AD + Qt + c 




■ X 


_ I2t if + 3t + c) 
I + D (4^3 Qt + c)^ 


(46) 


which thus explicitly gives the trajectories in the x —Om-plane in parametrized form, 
where the constant c characterizes the diherent solutions, with c = 0 yielding the 
separatrix subset, i.e., the heteroclinic orbits from FL^ to PLi in S (a heteroclinic 
orbit is a solution trajectory that connects two distinct hxed points). The solution 
space structure is depicted in Fig. [see also the reduced state space in Fig. |l(b)j. 


Let us now consider the interior solutions in Fig. 4(b) with A = \f3j2 in an him — q 
diagram and in several cosmographic diagrams (see Fig.[^, together with the ACDM 
model with = 0.3 and Dao = 0.7. To have a feeling for the evolution in terms 
of redshift one can compare with the graphs in Fig. which gives a good estimate 
for the history also for the other solutions when expressed in the redshift z, as long 
as they are close to the ACDM model. As can be seen, these models are deviating 
rather signihcantly from the ACDM model, and therefore a considerably smaller 
value of A is needed to obtain an evolution that is close to ACDM dynamics, even 
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(a) The scalar field boundary flm = 0; A = (b) The state space for dust and A = 

Figure 4: State space structures for dust and V = Voexp(—A0) with A = \f3fl. 


when Vtm is close to 1 initially. Furthermore, in a certain sense, discussed next, the 
present class of models is as special as the scalar held models with a cosmological 
constant, and they are therefore unlikely candidates for solving any of the issues one 
might have with ACDM cosmology. 


3.3 Symmetries of frozen ACDM dynamics 

It is only the cases of a constant or an exponential potential that effectively reduce 
the problem from three to two dimensions by decoupling the equation for the scalar 
held variable. The underlying reason for this is that these models admit scaling 
symmetries. 

In the case of an exponential potential, a translation of the scalar held leads to 
a scaling of the potential, which in combination with an appropriate scaling of t 
yields a scaling of furthermore, a scaling of the spatial coordinates leads to a 
scaling of the scale factor a and thereby also of pm in the case of a linear equation of 
state with ^ 0, which results in a scaling symmetry of Einstein’s held equations. 
This in turn gives rise to a one-parameter set of equivalent solutions, being the 
reason for the decoupling property, and it is also the underlying reason for why 
there exist “scaling solutions,” which more appropriately should be referred to as 
homothetic self-similar solutions since the spacetime geometry of these solutions 
admits a homothetic Killing vector held. For models with both a perfect huid with 
a linear equation of state and a scalar held with an exponential potential, a scaling 
solution can only exist if pm and behave in the same manner, which requires 
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Figure 5: Diagrams for the observable Qm and the cosmographic parameters H, 
q, j (the Hubble, deceleration, and jerk parameters) for dust and a pot ential V = 
Vq exp(— -^3/2 0) for the solutions (with the same notation) in Fig. 4(b) The thick 
line represents ACDM cosmology with = 0.3, Daq = 0.7. The horizontal lines 
with flrn = 0.3, H/Hq = 1, and q = —0.55 represent the present situation, i.e., 
z = 0. 


7 m = 7 (/>- The existence of a hxed point EM on the reduced state space x — Ftm 
therefore necessarily depends on this feature, but it turns out to not always be 
possible since EM only exists when A > ydF)^; furthermore, the global attractive 
property of EM, which does not follow from a hxed point analysis, follows from a 
monotone function that can be derived from another kind of symmetry of the held 
equations that is associated with coordinate scalings (this is a general mechanism, 
discussed for anisotropic models and perfect huids in IZIl)- For a constant potential 
scalings of the coordinates instead bring a solution with a given V = A to a solution 
with a diherent A, and because A carries dimension (length”^), there are no solutions 
admitting a homothetic symmetry when A 7 ^ 0; i.e., this scaling symmetry is not 
quite the same as that in the exponential case. 

The existence of symmetries is what makes these models tractable, and therefore 
also popular. However, the symmetries endow these models with special properties, 
which suggests that perhaps it is not a good idea to consider these models as role 
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models, since other models do not admit similar symmetries and properties, or at 
least that some carefnl considerations may be needed. Moreover, originally, before 
one realized that the deceleration parameter q is evolving, it was thonght that scale 
invariance and thereby a constant q was qnite desirable, since this wonld solve 
hne-tnning problems snch as coincidence and energy scale problems. On the other 
hand, the symmetries also make the models qnite inflexible. Dne to the changing 
evolutionary history of q, where q is currently negative, the scaling symmetries must 
be broken, while still leading to models that in some way at least alleviate various 
hne-tnning problems; for these reasons, and from a phenomenological perspective, 
scale-invariant models are no longer the favorite candidates as dark energy models. 


4 Dynamical ACDM dynamics 


4.1 Inverse power-law potentials 


In the inhuential papers pa Peebles and Ratra argued that dark energy can be 
phenomenologically modeled by a minimally coupled scalar held with a potential 
with a shallow tail, which results in the dark held energy density decreasing more 
slowly than the matter energy density to its “natural” value - zero. In particular 
they considered an inverse power-law potential, 

V = 0>O, Q!>0 


In addition they assumed a matter dominated hat FLRW universe after inhation, 
which resulted in a set of approximate hat FLRW equations for which they found a 
particular solution, which was shown to be linearly stable within the approximate 
context. The associated solution of the exact FLRW equations has subsequently 
been referred to as an attractor or tracker solution [5] , and has resulted in numerous 
papers. In [12] a thorough local dynamical systems analysis of the present models 
was performed, based on previous work, e.g., |22l |23l El] (for additional references 
see, e.g., mm)- Here we will use a slight variation of this approach and use the 
following bounded scalar held variable: 


Z 


1 

1 + A 


0 

(j) + a 


A 


1 - Z 

z 



(48) 


where Z thereby is monotonically increasing in 0 (in [T3] 1 — Z was used as the 
scalar held variable). 


The present case has an unbounded A, and therefore we use the dynamical sys¬ 
tem (13), where we choose g{Z) = Z, and thus gx_ = 1 and 




dr 

df 


dt 

df 


H-^Z. 


(49) 
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and the following dynamical system for the state vector {x,Qrn, Z)\ 

— = -(2 - q)xZ + ^ -{1- Z) (l- x^ - Dm) , (50a) 

^ = 3 [2a;2 - (1 - 12^)] D^Z, (50b) 

— = (1 — Z)'^ Zx, (50c) 

dr a 

where 2 — q = 3(1 — x^) — |f2m- 

Dne to the discnssion abont global dynamics in Sec.[^ it follows that the conclnsions 
obtained from the local hxed point resnlts in |T3] correspond to global asymptotic 
featnres. Hence a two-parameter set of solntions originates from M^, and a one- 
parameter set from FL^, where the tracker solntion originates from FLq. Toward the 
fntnre all orbits on S and end at the global fntnre attractor dSij^ The solntion 
strnctnre for a = 6, and especially that of the snbset that originates from FL^, is 
depicted in Fig. 



(a) The {x,il.rn,Z) state space. (b) Projection of state space onto 

(x, Hm)- 

Fignre 6: Solntion strnctnre for the state space for models with dnst and a scalar 
held with inverse power-law potential V = with a = 6. 

®Since lim 0 _).+oo F((^) = 0 it follows that the future asymptotic state dSi now corresponds 
to the Minkowski spacetime instead of the de Sitter spacetime. However, because the matter 
energy density goes to zero faster than the scalar field energy density it still follows that g —>■ — 1 ; 
the asymptotic evolution thus resembles that of the de Sitter spacetime, since asymptotically the 
spacetime approaches that of a Minkowski spacetime in a (local) foliation where g —>■ — 1. A similar 
statement holds for the models in the next subsection. 
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It is of some interest to complement the dynamical systems picture with a heuristic 
description where we regard the scalar field as a particle moving in a potential, while 
losing energy due to the friction force —3H(j). Due to the existence of a potential 
wall that is so steep that A —)■ oo when 0 —)■ 0, all “scalar field particles” (solutions) 
either come from 0 —>■ — cx) and bounce against the potential wall (the two-parameter 
set of solutions from M]“) or from being initially still at some 0 (no initial kinetic 
scalar held energy, i.e., a; = 0) and then rolling down the potential in the positive 
0-direction; the tracker solution corresponds to the solution that starts by being 
initially still with an inhnite scalar held energy associated with the limit 0 —0 
(note that this result also follows from the asymptotic explicit results concerning 
this solution given originally in [SI 0]). 

All solutions therefore eventually move in an increasingly shallow potential, slowed 
by “friction,” and hence with decreasing scalar held energy, although the matter 
energy decreases faster, thus leading to a future asymptotic (quasi) de Sitter state 
at 0 —>■ d-cx) (i.e., dSi). 

Loosely speaking, the tracker solution is connected with a bifurcation associated 
with going from a potential with hnite A to inhnite A, transforming the attracting 
focus of the exponential case with large A to a center. The orbits on Z = 0 are 
characterized by Dm = const. They therefore correspond to the periodic orbits 
found in the monomial case at late times, given in diherent variables in [T3], but in 
the present case they are interrupted by being cut in half by the line of hxed points. 
Mg, which is a formal consequence of the new time variable. The physical reason 
for this feature is that in contrast to the monomial case there are no scalar held 
oscillations at late times for inverse power-law potentials. 

As can be seen from Fig. the tracker solution forms part of the boundary of the 
FL^ saddle subset, which therefore is not a separatrix surface. Because of its initial 
linear stability property in the Peebles-Ratra formulation [310], which was seen as 
alleviating the coincidence problem, nearby solutions are attracted to it. However, 
there exists an open set of solutions that are further away from it that are not, 
even if Qm is close to 1 initially (those with fairly large initial Z). This is to be 
expected since the Peebles-Ratra stability analysis is only local. It should also be 
pointed out that there is no a priori reason for believing that the local attracting 
property holds globally in timej^ However, all solutions near the tracker solution, 
indeed all solutions that are near or on the subset that originates from FL^, are 
eventually shadowing the center manifold of the hxed point dSi on and it is 
the center manifold solution on that is an “attractor solution” rather than the 
tracker solution, since any solution close to this solution is attracted to it because 
it is a center manifold, but note that it is dSi that is the formal future attractor. 

®See, e.g. m for an example of an attracting solution that is locally but not globally attracting. 

^Incidentally, it is the center manifold solution on that the slow-roll approximation, given 
in, e.g., 0, describes approximately; see [Tsi [n] for a similar situation at early times instead of 
late times. 
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As seen from Fig. the tracker solution and the solutions that are near to it on 
the subset that originates from FL^ are those solutions on this subset that diverge 
from ACDM dynamics the most (in this sense they are the observationally worst 
solutions in this subset of solutions, which is also seen by plotting the solutions in 
diagrams for the observables; see Fig. and quite substantially when a = 6. Let 
us therefore consider the dynamics for models with lower a. 



(a) Urn — <7 diagram. 



(c) h — j diagram. 





(b) h — q diagram. 


/ 





0 0.2 0.4 0.6 0.8 1 

J 

(d) q — j diagram. 


Figure 7: Observables for dust and an inverse power-law potential V = Vo0“" 
with a = 6. The thick solid line represents ACDM cosmology; the horizontal lines 
represent now, i.e., zero redshift; and the other lines correspond to the solutions in 
Fig.§ where the grey line corresponds to the tracker solution. 

The solution structure for a = 1/2, and especially that of the subset that originates 
from FL^, is depicted in Fig. The topological structure is the same as for a = 6, 
but the subset that originates from FL^ is now located closer to the x = 0-plane 
and as a consequence the dynamics of solutions with Qrn close to 1 initially more 
closely resembles that of the ACDM models. These results are therefore in line with 
the conclusion in [6], that it is only for quite small a that the tracker solution is 
observationally viable, which restricts the observational relevance of these models. 
The deviations from ACDM cosmology can be explicitly seen in the diagrams for 
the observables in Fig. 

Note that Figs. and also show that it is not only the tracker solution that attracts 
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(a) {x,V,jn,Z) state space. (b) Projection of {x,ilm,Z) state space onto 

(x, 

Figure 8: The solution space of dust and an inverse power-law potential V = 
with a = 1/2. 


nearby solutions, but the whole FL^ saddle subset. The underlying reasons for this 
are that (i) the whole line of fixed points FL^ acts as a line of saddle points, and (ii) 
all solutions in this subset are attracted to the center manifold of dSi at late times 
before asymptotically approaching the global future attractor dSi. 

Finally, it is worth noting that the underlying reason for why the early behavior of 
the tracker solution exhibits a “tracking property” is because the present case can 
be viewed as a bifurcation where EMq merges with FLq in the limit A —)■ oo. The 
property that 7<^ = 7m = 1 for EM, irrespective of A (as long as EM exists until 
it merges with PL), can subsequently be exploited by comparing the reduced two- 
dimensional dynamics of a sequence of reduced exponential scalar field state spaces, 
as done in [131 ESI |23] , to yield an approximate heuristic description of the behavior 
of the tracking solution. Nevertheless, it should be clear that the dynamics of an 
inverse power-law potential model is globally very different from that of a given 
exponential potential, which is what can be expected from a model that does not 
exhibit scaling symmetries. 

The observational viability difficulties of the dust and inverse power-law potential 
models, including those of tracker solutions, arise because A —)■ cxd when </>—)■ cx), 
since a large A makes the dynamics deviate from ACDM dynamics. This motivates 
us to consider other types of scalar field potentials for which A is globally and asymp¬ 
totically regularized. For this reason, we therefore phenomenologically regularize A 
of the inverse power-law potential and use ([^ to derive V, in arguably the simplest 
possible way, so that A becomes finite and one obtains continuous deformations of 
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(b) h — q diagram 
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(d) q — j diagram 


Figure 9: Observables for dust and an inverse power-law potential V = with 

a = 1/2. The thick solid line represents ACDM cosmology; the horizontal lines 
represent now, i.e., zero redshift; and the other lines correspond to the solutions in 
Fig. 1^ where the grey line corresponds to the tracker solution. 


ACDM cosmology. 


4.2 A-regularization of inverse power-law potentials 


To obtain a type of potential that has a regular A and still behaves like an inverse 
power-law potential for large 0, we generalize the inverse power-law potential by 
choosing 




•^max 

T^+WW’ 


C>0. 


(51) 


Here Amax = a/C is the maximum value of A(0) while C describes the peak width 
of the present Using ([^, this “regularized” inverse power-law A(0) leads to 

the potential 


v = v„(^+^/c^ + 4>T“■ 


(52) 


It follows that lim,ji^±oo(</>^“U) = const. Note that this class of potentials contains 
the inverse power-law potential by setting C = 0 (restricting dynamics to </> > 0) 
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and a cosmological constant by setting a = 0. As an example, we give the potential 


and A for a = 6 and C = 10, depicted in Fig. 10 




(a) Potential diagram (b) A diagram 

Fignre 10: The potential V = V"o(0+ \/C‘^ + and associated X = a/+ cjY 
for a = 6 and C = 10. 


The above is (so far) a pnrely phenomenologically motivated potential, bnt so was 
the original motivation for the inverse power-law potential pa. The point here 
is to give a simple specific example of continnons ACDM scalar field deformations, 
illnstrating featnres that any observationally viable model arising from some fnnda- 
mental theory mnst have. 

To continne we choose 


Z = 


1 + ^’ 


Z = <j)+ VC'2 + 02, 


^2 _( J 2 ^2 _ 


2z 


2Z{1 - Z) 


(53) 


as the scalar field variable (here z shonld not be confnsed with the redshift). This 
leads to 

(54) 


z2 + c2(i - zy 


and the regnlar dynamical system 


x' = -(2 - q)x -b \l ^A(Z)(1 - 


Z' = y/Qa-^XZ{l -Z)x = VQ 


2Z\1- Zf 
Z2 + C2(1 - Z)2 


X 


(55a) 

(55b) 

(55c) 


where 2 — g = 3(1 — x^) — The qnalitative properties of this system follow from 
the general discnssion in Sec. We give an explicit representation of the solntion 


space for the case C* = 10 and a = 6 in Fig. The associated diagrams for physical 
observables are given in Fig. [T^ As can be seen the A-regnlarized models give mnch 
better resnlts than the corresponding models with an inverse power-law potential 
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(a) {x,V,jn,Z) state space. 


(b) Projection of state space onto 

(x, Urn)- 


Figure 11: Depiction of the solution space of flat FLRW models with dust and 
a scalar held with a potential V = Vo{(j) + and associated A = 

a/+ (fA for a = 6 and C = 10. 


with the same a, when C is chosen so that it is of the same order of magnitude as a, 
and even better if C is larger (cf. Figs. [^and[^ which both have a = 6). The results 
basically linearly improve with increasing C for a hxed a, which should not come as 
a surprise since Amax = a/C; for a = 6 and C = 100, and therefore Amax = 0.06, the 
dynamics for models with Qm close to 1 initially is almost indistinguishable from 
ACDM dynamics. Toward smaller C, the models increasingly yield similar results 
as the inverse power-law potential, which again is to be expected. 


4.3 Observational A conditions for continuous ACDM defor¬ 
mation 

The structure of the flat FLRW models with dust and inverse power-law potential 
strongly suggests that it does not suffice to just have a monotonically decreasing 
potential with a shallow potential tail, since even if A —)■ 0 when (j) —)■ -|-cx) a large 
X{(j)) for some (j) generally results in considerable deviations from ACDM cosmol¬ 
ogy. There are signihcant differences in the solution structure if A is globally and 
asymptotically bounded or not. 

This suggests that observational viability implies global and asymptotic boundedness 
conditions on A, so that the whole subset that originates from FL^ into S only 
deviates moderately from the x = 0 ACDM surface. It should be pointed out that 
this not only produces an observationally viable set of models, but also alleviates 
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(a) Hm — Q diagram 


(b) h — q diagram 




(c) h — j diagram 


(d) q — j diagram 


Figure 12: Observables for flat FLRW models with dust and a scalar field with the 
potential V = Vo(0 + \/C‘^ Y and associated A = a/\fCPM-~^ for a = 6 

and C = 10. The thick solid line represents ACDM cosmology; the horizontal lines 
represent now, i.e., zero redshift; and the other lines correspond to the solutions in 
Fig.[n] 


the coincidence and possibly other hne-tuning problems since all orbits with initial 
values of Qm close to 1 behave similarly because the invariant subset that originates 
from FL^ then constitutes an attracting separatrix surface, eliminating the need for 
any measure]^ 

^Because open sets of solutions do not follow the tracker solution closely, it was argued in m 
that solutions that were close to the tracker solution were in some sense generic due to the center 
manifold structures. However, such structures can be eliminated or changed by a change of vari¬ 
ables, while the fact that there are open sets of solutions that behave quite differently than the 
tracker solution in terms of physical observables cannot. It is therefore a question of measures if 
trackerlike behavior is generic or not, and a consensus about measures seems rather unlikely. On 
the other hand, if center manifold structures reflect behavior in physical observables, they might 
very well play an important role in cosmology. 
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A scalar field and two fluids 


In this work we have considered a matter content that consists of a fluid with¬ 
out pressure. However, the early Universe after inflation is radiation dominated and 
therefore somewhat more sophisticated models have both radiation and dust as mat¬ 
ter content. In a dynamical systems context this situation can be treated as follows. 
Gonsider two noninteracting fluids with energy densities pi and p 2 , respectively. It 
is then convenient to introduce 

pm = pi+p2, X=—, (56) 

Pm 

as the two perfect fluid variables instead of pi and p 2 f\ Assume further that the 
fluids obey linear equations of state characterized by the constant equation of state 
parameters 71 and 72 , where we without loss of generality set 72 > 71 (if 72 = 71 
we treat the problem as a single fluid). The system (|^ holds for fluids with general 
barotropic equations of state, which includes the present case. However, having two 
fluids instead of one with linear equations of state makes it convenient to augment 
the system (|^ with the additional variable y, which leads to 


x'=-{2 - q)x + (57a) 

= 3 [2x2 - 7,^(1 - G^)] (57b) 

Z' = Vq^x, (57c) 

a<p 

x' = 3(72-7i)x(1-x), (57d) 


where, as before, Gy = 1 — x‘^ — Ilm and flm = pm/^H^- The deceleration parameter 
is given by 


g = -1 + 3x2 -K = _1 + 3aj2 (^^^ + ^ 2(1 - x)) Hr, 


(58) 


since 7 ^ = = 71 X -f 72(1 — y). Due to the similarity with the single fluid 

case we choose to include the same boundaries, but with the addition of y we also 


®For alternative ways of handling several fluids, see, e.g., [251 [2S] . 
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consider x ^ [O 5 1]; where x = 0 (x = 1) corresponds to the single flnid case for p 2 
(Pi). 

Since 72 > 71 it follows that x is strictly monotonically increasing when 0 < x < 1- 
Fnrthermore, x can be solved in terms of the scale factor a by using the conservation 
equation dpi/dliaa = —S^iPi, which yields pi oc and therefore 


1 


1 

1 + A^exp[( 7 i - 72 )r|’ 


(59) 


where k = P 20 /P 10 = It follows that x 0 (x —t 1 ) when a —)■ 0 , i.e., 

r — )■ —(X) (a —)■ 00 , i.e. r —)■ +cx)). Thus x = 0 in the asymptotic past and x = 1 in 
the asymptotic future, which tells us that the fluid with the most soft (stiff) equation 
of state dominates asymptotically over the other fluid at late (early) times. As a 
consequence of the above features it follows that asymptotic future (past) dynamics 
reside on the x = 1 (x = 0) subset. 

Because it is possible to solve the conservation equations for each of the two fluids in 
terms of the scale factor, this gives rise to a constant of motion, which also reflects 
that the problem is really a three-dimensional one. To obtain this conserved quantity 
we hrst note that 

rii = xfim, ^2 = (1 - x)I^m- (60) 

Then, since pi oc and p 2 oc it follows that 


72 QT'2 

const. = % oc ^(1/2)72-71 = ^72(1 _ ;^)-7lf]72-7i(^2)72-7l^ (Ql) 

P 2 ii2 


where is given in terms of the state space variables in Eq. (19). 


The matter subset Dm = 1, i.e., a: = 0, Dy = 0, now corresponds to the two- 
dimensional system 


Z' = 0, (62a) 

x' = 3(72-7i)x(1-x), (62b) 

which has the solution Z = const., x = 1/(1 + ^exp[( 7 i — 72 )t]). The analysis of 
the single fluid case is directly translatable to the present context, and the line of 
hxed points FL^ is now replaced with the above two-dimensional subset. It follows 
that the inclusion of radiation changes little of the previous discussion as regards 
qualitative dynamical features. 

It should perhaps be pointed out that instead of x and Dm it is possible to use 
f2i = pi/ 3 if 2 as variables, where 

7 ^ 3 

—— = [2g — ( 37 * — 2)]flj, q = —1 -|- 3x^ + -( 7 if 2 i -|- 72122 ), (63) 

cLt 2 


and where fly = 1 — — f 2 i — f 22 . 
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